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DERIVATION OF FORMULAS FOR ROOT -MEAN-SQUARE ERRORS 
IN LOCATION, ORIENTATION, AND SHAPE IN TRIANGULATION 
SOLUTION OF AN ELONGATED OBJECT IN SPACE 

By Sheila Ann T, Long 
Langley Research Center 

SUMMARY 

Formulas are derived for the root- mean- square (rms) displacement, slope, and 
curvature errors in an azimuth- elevation image trace of an elongated object in space, 
as functions of the number and spacing of the input data points and the rms elevation 
error in the individual input data points from a single observation station. Also, formulas 
are derived for the total rms displacement, slope, and curvature error vectors in the tri- 
angulation solution of an elongated object in space due to the rms displacement, slope, 
and curvature errors, respectively, in the azimuth -elevation image traces from different 
observation stations. 


INTRODUCTION 

The purpose of this paper is to relate the errors in location, orientation, and shape 
in the triangulation solution of an elongated object in space to the rms elevation errors in 
the individual input data points. An elongated object, in this paper, is defined as an entity 
whose lateral dimensions are so small in comparison to its length that it can be regarded 
as a curved line. One example of such an entity in space is an ionization trail resulting 
from an object entering the upper atmosphere. A second example is an ionized cloud 
released in space, such as the barium ion cloud released in the magnetosphere on 
September 21, 1971 (ref. 1). 

In reference 2 a formula is derived for the total distance (i.e., displacement) error 
vector, which specifies the error in location, in the triangulation solution uf a smoke trail 
in space. For curved lines in space, however, one also needs to specify the errors in 
orientation and shape. For these, formulas are needed for the total slope and curvature 
error vectors, respectively. In this paper the total rms slope and curvature, as well as 
displacement, error vectors are determined. 

This paper is composed of two main parts. In the first part a smoothed azimuth- 
elevation image trace through the azimuth-elevation input data points from a single 
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observation station is obtained by a least-squares procedure. (The input data points are 
taken from photographs of the elongated object in space and are initially in terms of right 
ascension and declination coordinates* A coordinate transformation is then made from 
right ascension and declination to azimuth and elevation.) Formulas are derived relating 
the rms displacement, slope, and curvature errors in this smoothed azimuth-elevation 
image trace to the rms elevation error in the individual input data points. 

In the second part of this paper, the errors in the triangulation solution of an elon- 
gated object in space are considered. These errors are manifestations of triangulating 
on azimuth-elevation image traces (and, consequently, input data points) from several 
observation stations that contain errors. If the input displacement, slope, and curvature 
errors in the azimuth-elevation image traces are assumed to be small, then one can 
assume that they give rise to the total rms displacement, slope, and curvature error vec- 
tors, respectively, in the triangulation solution. Formulas are derived relating the total 
rms displacement, slope, and curvature error vectors in the triangulation solution of an 
elongated object in space to the rms displacement, slope, and curvature errors, respec- 
tively, in the azimuth- elevation image traces. 

The total rms displacement, slope, and curvature error vectors are useful for 
determining the relative merits of two or more different triangulation procedures appli- 
cable to elongated objects in space. They are also useful for determining the effect of 
various experimental parameters — such as the number and relative location of the obser- 
vation stations, the spacing of the input data points along the azimuth-elevation image 
traces, etc. - on the triangulation solution of an elongated object in space. 

SYMBOLS 

coefficients of second-degree curve (n-1)^, n^, (n+ if 

random curvature, slope, and displacement errors in segment (n-1)^, 
n^, (n+1)^ 

d N> d W d pN ,d i^N total, east-west, radial, and north-south output displacement errors 

in solution curve at point N, due to error e^ 

dj/j^ displacement error in curve (n-l)^, n^, (n+t)^ 1 at = 0 

total output displacement error vector in solution curve at point N, due 
to error ej^ 


Ap gM Ql- 1 - 

a N’ c N’^N 


a P- b M 
a N’ D N’ C N 
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N,N^ 


n,n^ 

r N 

S,S P 

S,S P 

P 

P2 ,p 4 

tP t p 
iv N-I 


N ,l N 


u 


nm 


input displacement error in segment (n-l) p n^, (n+l)* 1 

distance between points n+1 and (n+l) p ' in input slope error 
determination 

distance between points n and n p ' in input curvature error 
determination 

coordinates along ordinate axis (i.e., elevation) of individual input data 
points 

unit vectors centered at point N in directions of increasing east longitude, 
geocentric radius, and geocentric latitude 

quantity which is directly proportional to number of input data points 

index for sequential labeling of individual input data points 

indices for sequential labeling of points along curves S and 

indices for sequential labeling of points along traces s and 

input radius of curvature error in segment (n-l)^, n p , (n+l)^ 

reference and triangulation solution curves of elongated object in space 

azimuth-elevation error -free and image traces from station p 

total number of input data points 

quantities depending on number of input data points 

> > 

unit vectors in directions of vectors N^N+l and N-ljN^ - 

> > 

unit vectors in directions of vectors N,N+1 and N,(N+1)^ 

coordinates along abscissa axis [I.e., azimuth x cos(elevation)] 
of individual input data points 
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v^w* 1 continuous variables along abscissa (i.e., azimuth x cos(elevation)] and 

ordinate (i.e., elevation) axes of second-degree curve (n-1)^, n^ 
(n+l) M 


rectangular coordinate system with origin at center of earth, x axis in 
equatorial plane toward Greenwich meridian, z axis toward north, 
and y axis to form right-hand orthogonal triad 


4i'^lN >,3 2N 
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M 

nm 


c 


n 

N 


V 


N 


0 N ,0 1N’ 0 2N 


input slope error in segment (n-1)^, n^, (n+l) fI 

total, parallel, and perpendicular slope errors in solution curve at point N, 
due to error 

spacing between points n and n+1 

spacing between coordinates u^ m 

deviations associated with individual input data points 

vertical displacement of curve (n-1)^, n^, (n+1)^ at vj^ = 0 

curvature of cui've (n-1)^", n^, (n+1)^ at vff = 0 

total, parallel, and perpendicular rms curvature errors in solution curve 
at point N 


^N’^1N’^2N total, parallel, and perpendicular curvature errors, per unit input curva- 
ture error, in solution curve at point N 


i 


N 


*N’MN' k 2N 



input curvature error in segment (n-1) , n^, (n+1)^ 

total, parallel, and perpendicular output curvature errors in solution curve 
at point N, due to error 

total output curvature error vector in solution curve at point N, due to 
error 


A 


N’ a N 


east longitude of points N and N^ 1 
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M 


particular observation station 


■N’~AN’~pN’~</N 


CM eP- tM tP- 

? N ,? AN>Vn ,§ <£N 


total, east-west, radial, and north-south rms displacement errors 
in solution curve at point N 

total, east-west, radial, and north- south displacement errors, per 
unit input displacement error, in solution curve at point N 


P N’ P N 

°N 


geocentric radius of points N and N p 

rms elevation error in individual input data points 


°AN»°BN» °CN 


rms errors in A p , B p , and C p 


°aN»°bN ,cl cN rms curvature > slope, and displacement errors in segment (n-l) p , 

n p , (n+1 ) p 


N 


# 


geocentric latitude of points N and N p 
slope of curve (n-l) p , n p , (n+l) p at v p = 0 


^N , '^lN ,fi 2N total > parallel, and perpendicular rms slope errors in solution curve at 
point N 

• W N» W 1N ,C0 2N parallel, and perpendicular slope errors, per unit input slope error, 

in solution curve at point N 


ANALYSIS 


Relationship Between rms Displacement, Slope, and Curvature 
Errors in an Azimuth-Elevation Image Trace From a 
Single Observation Station and rms Elevation 
Error in Individual Input Data Points 

In figure 1 an azimuth -elevation plot from the observation station \i of an image 
trace s p and an error-free trace s is shown. An azimuth-elevation image trace 
approximates the actual photographic image of an elongated object in space. An azimuth- 
elevation error-free trace is the trace that would result if no errors whatsoever occurred 
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from the data acquisition system and the data reduction procedure. The image trace seg- 
ment (n-1)^, n^, (n+1)^ differs from the error-free trace segment n-1, n, n+1 in 
location (i.e., displacement), orientation (i.e., slope) as defined by the first-order deriva- 
tive, and shape as defined by the second and higher order derivatives. In this paper the 
highest order derivative considered will be the second; hence, shape will be defined by the 
single parameter curvature. 

For convenience the azimuth-elevation input data points within the segment are 
selected so that they are equally-spaced along the [azimuth x cos(elevation)] coordinate 
axis. The segments are sufficiently small so that they can be accurately described by a 
second-degree curve. Each second-degree curve is obtained by a least-squares procedure. 
The segment (n-l)^, n^, (n+1 f', centered about the point n^, of the azimuth- elevation 
image trace s ^ is one such resulting second-degree curve. 

The coordinates u^ m [i.e., azimuth x cos(elevation)) and h^ m (i.e., elevation), 
where m = -j, -j+1, ..., -1, 0, 1, ..., j-1, j, of the individual input data points in the interval 
which gives rise to the second-degree curve (n-l) M , n M , (n+1)^ are shown in figure 2. 
The deviations associated with these individual input data points are 




B&u^ 

N nm 



-h^ 

nm 


(V) 


The coefficients A{J, bJJ, and are to be determined. The sum of the squares of 

the deviations is 


u 


m=-] 


/ ,, \ 2 i r 

/ > 

,2 I 2 


C = l 

x 7 m=-i ^ 


) + B fKm "CjJ - hfj 

(2) 


The following is a generalization of the procedure set forth in reference 3 for the 
best fit of a straight line to experimental data. In this analysis the best fit of a second- 

j / \2 

degree curve to the data of interest is the one which minimizes V fe^ m j . Taking the 

m=-j ' ' 

first partial derivatives of equation (2) with respect to A^, B^, and C^, respectively, 
and then setting each of these three resulting expressions equal to zero, one arrives at 



J 

R M V 
U . 
m=-j 





( 3 ) 
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a n E.« m ) + b S5 E.(«L) +C S 2.«L- = ° 


m=-J 


m=-j 


A £ Z.{<J + B S 

rvk— v ' r>"i — _ 1 


m E «*+»<*- E •*»-« 


/J u nm ^ i'-N Li . nm 
m=-j nv=-j 


The quantity p is the total number of input data points in the interval of interest. 


p = 2j + 1 (6) 

Without loss of generality the 2j + 1 equally-spaced coordinates u£ m can be centered 
on a local coordinate origin ^i.e., uj^, = 0^. Then, 


y = o 

'-i . nm 
m=-j 


1 / \2 
E , < m - »2 6 ‘ 

n=-r 

j / \3 

E . Km - « 

*=-J 

E.Km)W 

'VI— ' 


The quantity 6 is the spacing between the coordinates u£ m of the individual input data 
points. The quantities pg and p^, which depend on the number of input data points in 


the interval of interest, are 


= 2(l 2 + 2 2 + 3 2 + . . . + j 2 ) 


_ id + l)(2j + l) 
3 
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p 4 = 2(l 4 + 2 4 + 3 4 + . . . + j 4 ) 

_ .1 (j + l)(2j + l)(3,j 2 + 3j - l) 
15 


( 12 ) 


Substituting equations (7) to (10) into equations (3) to (5) where appropriate, one finds that 
equations (3) to (5), respectively, become 


£ . h nm( u nm) 


(13) 


m=-j 


P2 5 X- 


m=~] 


nm 


(14) 


P 2 6 2 A£ 


N + PCft' E 


m =~J 


nm 


(15) 


after rearranging. Solving the three simultaneous equations (13), (14), and (15), one 
obtains the coefficients A^j-, B^., and Cj^., 

A „ i-j[ P Kn.) 2 -P2^^m 


N 


(pp 4 - p!) 


(16) 


r M _ m=-j 

b n 


E ujf. 


nm nm 


5 2 p 2 


(17) 


r fi _ m=-j 


g.. [^4® 2 - P2(<m) 2 


hP 

nm 


: (pp 4 -pl) 


(18) 


One can assume that the individual input data points are uncorrelated. If the rms 
error in the individual coordinates h^^ is denoted by o^, then the rms error °AN 
in a n is obtained from the following: 
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/u x 2 jiKJ-v 2 ] W 2 

W = 1 —r; ^5 


using equation (16) 


(pp 4 - P^) 

_ (°n) [p 2 ( P 4 54 ) ~ 2pp 2 5 2 (p 2 5 2 ) +P2^ 4 (P)] 
5 8 (pp 4 - Pg) 2 


using equations (10), (8), and (6) 

2 


6 4 (PP 4 - P|) 

Hence, the rms error o^j in A^- is 


J 1 

M _ °N / P 
°AN = ^2' 


1/2 


(19) 


PP 4 - Pg 


Similarly, the rms errors in and Cj^, denoted by o^ N and o^ N , respectively, 


'BN CN’ 


are 


« 1 

°BN = 6 . 1/2 


(«•*) 


°CN ~ °n| 


P 4 


,1/2 


(20) 


( 21 ) 


PP 4 - 


From equations (19), (20), and (21), using equations (6), (11), and (12), one notes 
that, for large j, 


a AN varies approximately as — 


9 



1 


°BN varies approximately as 

•6* varies approximately as 

The quantity 6 is the spacing between the coordinates u^ m of the individual input data 
points, and the quantity 3 is directly proportional to the number of input data points in 
the interval of interest. 

The equation for the second-degree curve (n- 1 )^, n^, (n+l) fJ ‘, centered about the 

point n^, is 


1 — 

‘"V' 1 


\~ J / 


<> 1/2 


(i) 







(22) 


where v^ and shown in figure 2, are the continuous variables along the abscissa 
(i.e., azimuth x cos(elevation)] and ordinate (i.e., elevation) axes, respectively. The ver 
tical displacement ^ of this curve, at = 0 (i.e., at the point nM-) r is 




w, 


n 


Ivj^O 



The vertical displacement error d£j^ is 


- dC £ 

The normal displacement error di^ in the curve (n- 1 )^, n^, (n+ 1 )^ 1 , at vJJ = 0, is 

dw^ N 


- < cos I tan 


•1 n 




v^O 

n 


= dC^ cos^tan - ^ 




r / \ 21 

1/2 

L 1 + K)J 
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Therefore, the rms displacement error o^- in the azimuth- elevation image trace 
segment (n-l)^, n^", (n+1)^, as related to the rms elevation error in the individual 
input data points in the interval of interest, is 



where oj^j is given by equation (21). 

The slope i of the second-degree curve (n-1)^", n^, (n+1)^, at vj^ = 0, is 


n 




-tan-1 b M 
The slope error di/^j. is 



Therefore, the rms slope error crj^ in the azimuth-elevation image trace segment 
(n-l)^, n^, (n+l)^, as related to the rms elevation error in the individual input 
data points in the interval of interest, is 



where is given by equation (20). 

The curvature of the second-degree curve (n-1)^, n^, (n+1)^, at v^ = 0, 
is 
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Therefore, the rms curvature error in the azimuth-elevation image trace segment 

n^ 1 , (n+1)^ 1 , as related to the rms elevation error in the individual input data 

points in the interval of interest, is 



where and o^ N are given by equations (19) and (20), respectively. 


(25) 
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Relationship Between rms Displacement, Slope, and Curvature Errors 
in the Triangulation Solution of an Elongated Object in Space and 
rms Displacement, Slope, and Curvature Errors, Respectively, 
in Azimuth-Elevation Image Traces 

Separation of errors .- The output displacement, slope, and curvature errors in the 
triangulation solution of an elongated object in space are assumed to be small and linearly 
related to the input displacement, slope, and curvature errors, respectively, in the 
azimuth -elevation image traces. A linear error analysis model will be derived using 
perturbations of triangulation solutions to establish the required linear relationships. 

The problem reduces to defining the deviation of one curved line in space from a 
second curved line in its immediate neighborhood. The first curved line is called the ref- 
erence curve and is denoted by S. The reference curve S, shown in figure 3, is assumed 
to have the same position and form as the elongated object in space of interest. The 
points along the curve S are sequentially labeled ..., N-l, N, N+l, ... . The geocentric 
latitude, east longitude, and geocentric radius of the point N on the curve S are 
denoted by <^, X N , and p N , respectively. In computations one should space the points 
along the curve S at small, equal intervals in geocentric latitude, east longitude, or 
geocentric radius, depending on the position and form of the elongated object in space of 
interest. Also, the points along the curve S should be spaced so that the distances 
between the successive points are small compared to the geocentric radii of the points. 
Hence, the arc segments between the successive points along the curve S can be 
regarded as straight-line segments. 

From the curve S, one can derive the azimuth-elevation error-free trace from 
each of the observation stations of interest. The azimuth- elevation error-free trace s 
from the observation station p is shown in figure 4. The points n-l, n, and n+l 
along the trace s correspond to the points N-l, N, and N+l, respectively, along the 
curve S. 

The input displacement, slope, and curvature errors are assumed to be small. 

Hence, one can assume that they give rise to the total rms displacement, slope, and curva- 
ture error vectors, respectively. In other words, the total rms displacement, slope, and 
curvature error vectors are assumed to be effectively decoupled and, consequently, can 
be examined separately. 

One notes that configurations of one, two, and three points define zero, first, and 
second order derivatives, respectively, if the derivatives are expressed in finite difference 
notation. This fact will be used in this part of the analysis. In other words, configurations 
of one, two, and three points will be used in the determinations of the total rms displace- 
ment, slope, and curvature error vectors, respectively. 
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The general procedure for determining the total rms displacement error vector, 
for instance, will be as follows: Introduce an input displacement error into the error- 
free trace s from station p to form an image trace sM- ; triangulate , using the 
trace s M from station p and the error-free traces from the remaining stations, to 
obtain a triangulation solution curve S^; compare the curve Sl J - to the curve S to find 
the total output displacement error vector due to the input displacement error introduced 
at the station p; take the ratio of the output to the input displacement errors; and com- 
pute the rms sum of this ratio over all of the stations. Analogous procedures will also 
be carried out with regard to the total rms slope and curvature error vectors. 


Displacement errors .- The determination of the total rms displacement error vector 
will be considered first. For this the points n-1, n, and n+1 on the azimuth- elevation 
error-free trace s from the observation station p are displaced through radians 
along the normals to the trace s at the points n-1, n, and n+1, respectively, to form 
the points (n-1)^, n^, and (n+1)^, respectively, as shown in figure 4. The segment 
(n-1)^, n^, (n+1)^, extended, is the azimuth-elevation image trace from the obser- 
vation station p. The quantity is the input displacement error in the azimuth- 

elevation image trace segment (n-l)^, n^, (n+1)^ from the observation station p. 

The triangulation solution of the elongated object in space is then calculated, using 
for input data the points along the azimuth-elevation image trace from the observa- 
tion station p in conjunction with the points along the azimuth- elevation error-free 
traces from all of the remaining observation stations. Three appropriate triangulation 
procedures for a curved line in space are described in references 4, 5, and 6. 


The resulting triangulation solution curve of the elongated object in space is 
denoted by S 11 and is shown in figure 3. The points along the curve are spaced 
according to the same scheme as that used for spacing the points along the curve S. 

They are spaced at the same small, equal intervals in geocentric latitude, east longitude, 
or geocentric radius that were used for the spacing of the points along the curve S. The 
point on the curve which has the same geocentric latitude, east longitude, or geocen- 
tric radius as the point N on the curve S is denoted by The geocentric latitude, 

east longitude, and geocentric radius of the point are denoted by j-, A^, and 

Pj^, respectively. (One should notice that no correspondence exists between the points 

and n^, unlike the situation for the points N and n.) 
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The vector N,N+1, from the point N to the point N+l, shown in figure 3, is 


N ’ N+1 = P N+1 C0S %+l sin ( X N+l * X n)*X 

+ p N+l|®* n ^+1 cos ^ ~ COS ^N+l sp5 C0S ( X N+l " X n) 
P N+lj^ 0S %+l C0S % C0S ( X N+1 “ X n) + Sjjl %+l Sin % 


' P N V 


( 26 ) 


where i^, i^, and i are the unit vectors centered at the point N in the directions of 
increasing east longitude, geocentric latitude, and geocentric radius, respectively, where 
(^j, X^, and p^ are the geocentric latitude, east longitude, and geocentric radius, 

respectively, of the point N, and where <^ + ^, X N+ ^, and P N ^ are the g eocentric lat- 

itude, east longitude, and geocentric radius, respectively, of the point N+l. The unit 

a ^ 

vector t^j, shown in figure 3, which is in the same direction as the vector N,N+1, is 

> 


; _ N,N+1 

r N~ * 

N,N+1 


The vector N,N P , from the point N to the point N p , in figure 3, is 


(27) 


N,N P - cos 41 sin ^X^- - X^i-^ 

+ PjjJsin 41 cos - cos 41 sin ^ cos (x p - X N ) 
r 

+ < Pjl^ os 41 cos C0S ( X N ■ X n) + sin 41 sin 




P N > 1- 


(28) 


where 41> 4j, and p^ are the geocentric latitude, east longitude, and geocentric 

radius, respectively, of the point N p . 

The vector d^, shown in figure 3, which is drawn from and normal to the curve S, 
at the point N, to the curve S p , very closely approximates the shortest vector distance 

from the curve S to the point N p . The vector d^ is the total output displacement 
error vector in the triangulation solution curve at the point N, due to the input displace- 
ment error ej^. 
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From figure 3 one sees that 


dft = N,N^ - l N,I# • t Jt 


N 


(29) 


Substitution of equations (27) and (28) into equation (29) leads to an equation for the vector 
in the following form: 


3 N 


,11 * 

~ d AN l \ 


+ d $N k 


+ d pY 


(30) 


where d(^, d 0N’ and d pN are the east_west (i- e *> ^ the direction of increasing east 
longitude), north-south (i.e., in the direction of increasing geocentric latitude), and radial 
(i.e., in the direction of increasing geocentric radius) output displacement errors in the 
solution curve at the point N, due to the input displacement error e^. The total output 
displacement error d^ in the solution curve at the point N, due to the input displace- 
ment error e^, is 



Hence, the east -west, north-south, radial, and total displacement errors, denoted by 
^XN’ ^0N* ^pN» and ^N’ res P ec tively, per unit input displacement error, in the solution 
curve at the point N are 


_ a XN 
9 *N" u 
e N 


(32) 


_ q 0N 

p. 

e N 


(33) 





(34) 






( 35 ) 
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Now, the east-west displacement error for instance, in the solution curve at 

the point N is due only to the unit input displacement error in the azimuth- elevation 
image trace segment (n-1)^, n^, (n+1)^ from the observation station p. If a random 
normal displacement error, c^-, with zero mean exists in the segment (n-1)^, n^, 
(n+1)^" from the station p, then the east-west displacement error in the solution curve 
at the point N, due to the random displacement error cj^, is c n)* ^ ence > Ihe 

east-west rms displacement error E^ N in the solution curve at the point N, due to the 
random displacement error in the azimuth-elevation image trace segment (n-1)^, 

n^, (n+1)^ from each of the respective observation stations, is 


.-, 1/2 



The bar denotes the mean value, and the summation is over the total number of observa- 
tion stations. Because the random displacement errors c^ at the various observation 
stations are uncorrelated, equation (36) becomes 


3 1/2 



If one assumes that the mean-square random displacement error m the azimuth- 

elevation image trace segment (n-l)t 1 , n^ 1 , (n+l)^ from each respective observation 

2 

station is equal to radians [the mean-square displacement error in the azimuth- 

elevation image trace segment (n-1)^, n^, (n+1)^ as given by equation (23)), then equa- 

tion (37) becomes 



Therefore, equation (38) is the formula for the east -west rms displacement error in the 
triangulation solution curve at the point N, as related to the rms displacement errors 
°cN in the azimuth- elevation image traces. Similarly, the north-south, radial, and 
total rms displacement errors, denoted by E ^ N , S pN , and S N , respectively, in the 
triangulation solution curve at the point N, as related to the rms displacement errors 
°cN in the azimuth-elevation image traces, are 
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1/2 





\ 2 / 


\ 2 | 

W 




( 39 ) 


(40) 


(41) 


Slope errors .- In this section the relationship between the total rms slope error 
vector in the triangulation solution of an elongated object in space and the rms slope 
errors in the azimuth- elevation image traces is considered. For this determination a 
triangulation solution of the elongated object in space is calculated, using for input data 
the points along the azimuth-elevation image trace segment n, (n+1)^, extended, 
denoted by s^, and shown in figure 5, from the observation station ju. in conjunction 
with the points along the azimuth-elevation error-free traces from all of the remaining 
observation stations. One notes that the segment n, (n+1)^ has one extremity on the 
trace s. This is permissible, without loss of generality, because the total rms slope 
error vector has been assumed to be effectively decoupled from the total rms displace- 
ment error vector. (The point (n+1)^ in fig. 5 is not necessarily the same as the 
point (n+1)^ in fig. 4.) The resulting triangulation solution curve S ,J " of the elongated 
object in space extends through the segment N, (N+l)^, as shown in figure 6. (There is 
no correspondence between the point (N+l) Ji in fig. 6 and the point (n+l)^ in fig. 5 and 
the point (N+l)^ in fig. 6 is not the same as the point (N+l)^ in fig. 3.) 

The angle shown in figure 5, between the segment n, n+1 and the segment 

n, (n+l/' t is the input slope error in the azimuth-elevation image trace segment (n-1)^, 
n^, (n+l)^ from the observation station ju. The angle shown in figure 6, between 

the segment N, N+1 and the segment N, (N+l)^ is the total output slope error in the 
triangulation solution curve at the point N, due to the input slope error a^. The first 
step in this section is to determine the equations for the total output slope error /3^ 
and the input slope error 
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> 

The vector from the point N to the point (N+l)^, in figure 6 , is N,(N+1)^. 



where ^N+l are S eocen t r i c latitude, east longitude, and geocen- 

tric radius, respectively, of the point (N+l)^. The unit vector th-, shown in figure 6 , 

which is in the same direction as the vector N,(N+1)^, is 


> 

111 _ N,(N+1) m 

N= l^l 


(43) 


Also from figure 6 one sees that the total output slope error in the solution 

curve at the point N, due to the input slope error q^, is 



sin 


-II 


t xF 
t N C N 


(44) 


where the unit vector t N is given by equation (27). Now, the unit vector that is normal 


t„ x i 


, which is seen from figure 6 . The local 


*N X l p\ 


vertical tangent plane, here, is the plane passing through the earth's center and tangent to 
the reference curve S at the point N. The total output slope error can be 

resolved into two components, the parallel output slope error (parallel to the local 

vertical tangent plane) and the perpendicular output slope error /3^ (perpendicular to 
the local vertical tangent plane). From figure 6 one observes that, for small slope 
errors, the perpendicular output slope error and the parallel output slope error 

0 ^ in the solution curve at the point N, due to the input slope error are 
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nl 1 - (t X t^ S ) * ‘p/ 

r 2N "" \ N X *N/ li ! j 






N p) 
1/2 


(45) ■ (4jn)^ 


( 46 ) 


From figure 5 one sees that the input slope error in the azimuth- elevation 

image trace segment (n-1)^, n^, (n+1)^ from the observation station ju is 


= tan 


-1 


f b 

*N 


(As)ft 




N 


(47) 


The quantity (As)^ is the spacing in radians between the points n and n+1, and the 
quantity is the distance in radians between the points n+1 and (n+1)^. (in compu- 
tations one may set the quantity equal to the input displacement error ej^.) Hence, 
the parallel, perpendicular, and total slope errors, denoted by w^ N , and u>^, 

respectively, per unit input slope error, in the solution curve at the point N are 


0^ = ^ 
1N 

g n 

(48) 

e 

N>-p 

2! 

II 

(49) 

II 

¥ 

(50) 


Now, the parallel slope error u>^t, for instance, in the solution curve at the 
point N is due only to the unit input slope error in the azimuth-elevation image trace 
segment (n-1)^, n^ 1 , (n+1)^ from the observation station p. If a random slope 

error b^ with zero mean exists in the segment (n-1)^, n^, (n+l/ 1 from the sta- 
tion p, then the parallel slope error in the solution curve at the point N, due to the ran- 
dom slope error b^, is (^p^j b£|j.). The random slope errors at the various observation 
stations, like the random displacement errors, are uncorrelated. If one assumes that the 
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mean-square random slope error j^bj^ j in the azimuth-elevation image trace segment 

(n-1)^, n*\ (n+1)^ from each respective observation station is equal to (^n) radians 
(the mean-square slope error in the azimuth- elevation image trace segment (n-1)^, n^, 
(n+l) M as given by eq. (24)), then 


°1N “ 



(51) 


is the formula for the parallel rms slope error in the triangulation solution curve at the 
point N, as related to the rms slope errors in the azimuth-elevation image traces. 

Similarly, the perpendicular and total slope errors, denoted by and respec- 

tively, in the triangulation solution curve at the point N, as related to the rms slope 
errors in the azimuth- elevation image traces, are 




(53) 


Curvature errors .- In this section the relationship between the total rms curvature 
error vector in the triangulation solution of an elongated object in space and the rms cur- 
vature errors in the azimuth-elevation image traces is considered. For this determina- 
tion a triangulation solution of the elongated object in space is calculated, using for input 
data the points along the azimuth-elevation image trace segment n-1, n^, n+1, extended, 
denoted by s^, and shown in figure 7, from the observation station p in conjunction with 
the points along the azimuth-elevation error-free traces from all of the remaining obser- 
vation stations. One notes that the segment n-1, n^ - , n+1 has both of its extremities 
on the trace s. This is permissible, without loss of generality, because the total rms 
curvature error vector has been assumed to be effectively decoupled from both the total 
rms displacement and slope error vectors. (The point nM in fig. 7 is not necessarily the 
same as the point n^ in fig. 4.) The resulting triangulation solution curve of the 
elongated object in space extends through the segment N-1, N^, N+1, as shown in fig- 
ure 8. (No correspondence exists between the point N^ in fig. 8 and the point n^ in 
fig. 7, and the point N^ - in fig, 8 is not the same as the point N^" in fig. 3.) 
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The input curvature error •ft 1 »e., the reciprocal of the input radius of curvature 
error r^j in the azimuth-elevation image trace segment (n-l) w , n p , (n+l) p from the 
observation station ju is shown in figure 7. The total output curvature error vector 

in the triangulation solution curve at the point N, due to the input curvature error 

is shown in figure 8. The first step in this section is to determine the equations for 

the total output curvature error vector and the input curvature error 
> 

The vector N P ,N+1, from the point N p to the point N+l, in figure 8, is 


N P N+1 = 


r 


P N+1 C0S ^N+l Sm ( X N+l “ X n) ‘ P N C0S 4J sin ( X N " V)) 1 ; 


+ < | p N+l|^ ln ^+1 cos ’ cos ^N+l sm cos ( X N+1 “ X n) 
- pjjjsin cos ^ - cos sin cos(x{J - X N j 


P N+1 


sin ^ +1 sin + cos <^ +1 cos cos(x N+1 - X N ) 


P N 


sin sin (^ + cos cos cos (x p - X N ) 


(54) 


The unit vector T^, shown in figure 8, which is in the same direction as the vector 
N p , N+l, is 


_ N p N+l 

N "Fzr^t 

|n p ,n+i 
> 

The vector N-1,N P , from the point N-l to the point N p in figure 8, is 


(55) 
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N-l 


,nM = (^N C0S sin ( X N “ X n) + p N- 1 C0S ^N-l Sin ( X N " X N-l)) 1 X 
+ ^P^j^sin- <$• cos - cos ^ sin ^ cos (x p - \ N ) 

- P N _ijc°s 4>x sin - Sin cos cos (X N - * N _J^ ^ 

“ “I 

sin sin + cos ^ cos cos ^X.^j. - Aj^ 


+ < P N 


- P 


N-l 


sin sin + cos cos <^_ x cos(a n - A n _j) 


ip (56) 


where A^_p and are the geocentric latitude, east longitude, and geocen- 

tric radius, respectively, of the point N-l. The unit vector T^_j, shown in figure 8, 
which is in the same direction as the vector N-1,N P , is 


_ N-1,N P 

N - 1 = P^1 


If the curvature of the reference curve S is small compared to the curvature of 

the triangulation solution curve S p , then the total output curvature error vector in 
the solution curve at the point N, due to the input curvature error t is approximately 



(58) 


as is seen from figure 8. The total output curvature error Kjl in the solution curve at 
the point N, due to the input curvature error i^, is 



(59) 
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The total output curvature error can be resolved into two components, the parallel 
output curvature error (parallel to the local vertical tangent plane, as defined in 

the previous section on slope errors) and the perpendicular output curvature error 
(perpendicular to the local vertical tangent plane). The perpendicular output curvature 
error and the parallel output curvature error in the solution curve at the 

point N, due to the input curvature error are 



where the unit vector t^ is given by equation (27). 

Now, the point n ;i is approximately midway between the points n-1 and n+1 
along the azimuth-elevation image trace segment n-1, n^, n+1 and also the quantity 

g^, which is the distance in radians between the points n and n^, is much smaller than 
the input radius of curvature error r|^, as are shown in figure 7. Hence, one finds that, 
to sufficient accuracy, the input radius of curvature error rj^ in the azimuth -elevation 
image trace segment (n-1)^, n^, (n+1)^ 1 from the observation station ju is 



where (As)^j- is the spacing in radians between the points n and n+1. 

Hence, the input curvature error in the azimuth-elevation image trace seg- 
ment (n-1)^, n^, (n+1 )^ from the observation station p is 



(^In computations the quantity g^j- may be set equal to the input displacement error ej^.J 
Hence, the parallel, perpendicular, and total curvature errors, denoted by 0^, ^2N’ 
and 0^, respectively, per unit input curvature error, in the solution curve at the point N 
are 
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( 64 ) 



jJL 

eft = K Ji (66) 

N p 
l N 

Now, the parallel curvature error for instance, in the solution curve at the 

point N is due only to the unit input curvature error in the azimuth-elevation image 
trace segment (n-1)^, n^, (n+1)^ from the observation station p. If a random curva- 

ture error a^ with zero mean exists in the segment (n-1)^, n^, (n+1)^ from the sta- 
tion p, then the parallel curvature error in the solution curve at the point N, due to the 
random curvature error ajj, is ( b^ n ajj). The random curvature errors at the various 
observation stations, like both the random displacement and slope erro rs, ar e uncorre- 
lated. If one assumes that the mean- square random curvature error ^a£^ in the 
azimuth-elevation image trace segment (n-1)^, n^, (n+1)^ from each respective obser- 
vation station is equal to radians (the mean- square curvature error in the 

azimuth-elevation image trace segment (n-l)^, n^, (n+l) JX as given by eq. (25)), then 



is the formula for the parallel rms curvature error in the triangulation solution curve at 
the point N, as related to the rms curvature errors in the azimuth- elevation 

image traces. Similarly, the perpendicular and total rms curvature errors, denoted by 
© 2 ^ and © N , respectively, in the triangulation solution curve at the point N, as related 
to the rms curvature errors cr^ in the azimuth-elevation image traces, are 
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(69) 



CONCLUDING REMARKS 


Formulas have been derived for the rms displacement, slope, and curvature errors 
in an azimuth- elevation image trace of an elongated object in space, as related to the rms 
elevation error in the individual input data points. Also, formulas have been derived for 
the total rms displacement, slope, and curvature error vectors in the triangulation solu- 
tion of an elongated object in space, as related to the rms displacement, slope, and curva- 
ture errors, respectively, in the azimuth-elevation image traces. These total rms dis- 
placement, slope, and curvature error vectors specify the errors in location, orientation, 
and shape, respectively, in the triangulation solution of an elongated object in space. 
Therefore, the errors in location, orientation, and shape in the triangulation solution of 
an elongated object in space have been related to the rms elevation errors in the indi- 
vidual input data points. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., January 23, 1974. 
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Figure 1.- From observation station p, image trace segment (n-l)* 1 , n^, (n+1)^ 

differing from error-free trace segment n-1, n, n+l in location, orientation, 
and shape. 












